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ABSTRACT 


This  report  deals  with  the  basic  seismic  analysis  of  small  mag- 
nitude regional  Eurasian  events  observed  at  NORSAR.  The  analysis  was  severe- 
ly hampered  by  the  experimental  limitations.  However,  phase  energy  ratios 
(P/S)  appear  promising  as  a discriminant  despite  the  experimental  problems. 
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elusions  presented  are  those  of  the  authors  and  should  not  be  interpreted  as 
necessarily  representing  the  official  policies,  either  expressed  or  implied,  of 
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SECTION  I 


INTRODUCTION 


Teleseismic  signals  from  Eurasian  events  have  been  studied 
extensively  using  data  from  seismic  stations  and  arrays  throughout  the  world. 
Various  sophisticated  signal  enhancement  techniques  have  been  developed  to 
lower  detection  and  discrimination  thresholds.  However,  many  small  or  low 
magnitude  events  occur  which  are  not  detected  at  teleseismic  distances,  even 
with  the  present  processing  techniques.  Small  epicentral  distances  are 
a necessary  requirement  for  detection  of  these  events,  and  it  is  important  to 
have  knowledge  of  the  parameters  and  techniques  required  for  discriminating 
between  regional  small  shallow-earthquakes  and  explosions.  The  nature  of 
discrimination  at  short  distances  can  be  quite  d i'ferent  than  that  for  teleseis- 
mic distances,  because  of  the  differences  between  regional  and  teleseismic  m 's 

b 

and  because  of  the  great  variety  of  phases  available  at  regional  distances.  This 
report  considers  some  of  the  problems  in  regional  event  discrimination. 

In  tins  study  we  are  mainly  concerned  with  the  distance  range 
at  which  the  primary  wave  is  refracted  horizontally  along  the  Mohorovicic  dis- 
continuity (Pn).  Pasechnick  (1970)  indicated  that  the  distance  range  for  the 
first  arrival  (PJ  over  Eurasia  was  \ < 800  to  1200  kilometers,  which  he  terms 
'first-zone1.  The  distance  range  1200  £ A < 2000  kilometers  he  terms  the 
'second-zone'.  Carder  (1952)  on  the  other  hand,  uses  the  term  near  regional 
for  1501 A £650  kilometers  and  regional  650  < \<  1600  kilometers.  We  use 
the  terms  first-zone  and  regional  interchangeably  throughout  this  report  but  are 
always  referin6  to  distances  wilwin  the  general  range  of  150  1 4 < 2000  kilo- 
meters. 

First-zone  studies  in  Eurasia  have  been  reported  primarily  by 
Pasechnick  (1970).  He  observed  that  both  earthquakes  and  explosions  from  the 
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same  regions  had  the  same  phases  and  travel-times  but  that  dynamic  charac- 
teristics of  the  phases  were  different.  Laun  and  Becker  (1974)  investigated 
the  feasibility  of  using  NORSAR  data  to  examine  differences  in  the  dynamic 
characteristics  of  these  phases.  Their  overall  goals  were  based  on  results 
obtained  from  extensive  studies  by  many  authors  of  well  located  events  from 
the  Nevada  Test  Site  (NTS)  recorded  at  near  regional,  regional,  and  small 
telese’smic  distances.  However,  for  seismic  events  observed  at  NORSAR, 
precise  locations,  origin  times  and  depths  were  not  possible,  especially  for 
those  events  detected  and  located  by  NORSAR  alone.  Thus,  many  of  the  dis- 
criminant techniques  applied  to  NTS  events  were  not  successful  when  applied 
to  the  NORSAR  data.  This  does  not  mean  that  these  techniques  would  not  be 
effective  for  Eurasian  regional  events  observed  at  NORSAR  if  precise  hypo- 
center  parameters  were  known,  but  only  that  they  were  unable  to  apply  them 
with  any  degree  of  confidence  to  their  particular  data, 

Laun  and  Becker  did  find  one  promising  discriminant  despite 
their  experimental  problems,  and  that  was  the  P/S  energy  ratio.  Their  de- 
termination of  this  ratio  utilized  broad  time  windows  based  on  fixed  phase 
velocities  (Booker  and  Mitrinovas,  1964)  to  reduce  the  problems  due  to  inac- 
curate hypocentral  locations. 

In  this  report  we  have  enlarged  the  data  base  of  Laun  and  Becker. 
Unfortunately  no  additional  presumed  explosions  occurred  during  the  time  span 
of  this  new  event  ensemble.  However,  we  make  some  observations  and  suggest- 
ions concerning  the  experimental  problems  encountered. 

This  report  is  arranged  in  the  following  manner:  Section  II  pie- 
sents  the  data  base.  Section  III  deals  with  processing  and  analysis,  Section  IV 
presents  the  results,  Section  V contains  the  pertinent  conclusions  and  suggest- 
ions, and  Section  VI  gives  references. 
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SECTION  II 


BASIC  DATA 

The  NORSAR  array  at  Kjelier,  Norway  was  chosen  for  the 
source  of  first-zone  data  for  a number  of  reasons.  The  array  is  situated 
close  to  areas  where  presumed  explosions  have  occurred.  The  propagation 
path  for  these  signals  is  typically  continental  and  the  crustal  properties  in 
Scandinavia  have  been  studied  in  detail.  As  for  the  array  itself,  both  high 
quality  vertical  component  short-period  data  and  three  component  long-period 
data  are  recorded  there.  The  NORSAR  Data  Processing  Center  also  publishes 
a bulletin  of  events  detected  by  the  array.  This  bulletin  reports  many  low 
magnitude  regional  events  which  do  not  anpear  in  the  Preliminary  Determina- 
tion of  Epicenters  (PDE). 

These  NORSAR  data  were  obtained  in  two  ways.  The  short- 
period  data  are  ordered  from  the  Norwegian  Data  Processing  Center  at  Kjelier, 
Norway  and  mailed  to  the  Seismic  Data  Analysis  Center  (SDAC)  in  Alexandria, 
Virginia.  A time  lag  of  approximately  three  months  exists  between  the  data 
request  and  having  the  short-period  data  on  hand  for  processing.  All  long- 
period  data  are  transmitted  in  real  time  via  a trans- Atlantic  communication 
link  to  the  SDAC  and  recorded  there  on  magnetic  tape.  However,  data  quality 
problems,  combined  with  the  small  size  of  the  event  data  base  resulted  in 
having  only  five  events  with  usable  long -period  data.  Thus,  no  long-period 
data  results  are  presented  in  this  report. 

A first-zone  or  regional  event  list  was  compilrd  from  both  the 
PDE  bulletin  and  NORSAH  bulletin  for  the  period  January  1974  through  March 
1974.  Table  II- 1 lists  60  events  and  their  associated  event  parameters. 
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Designation 


TABLE  11-1 

EVENT  LIST  AND  PARAMETERS 
(PAGE  2 OF  3) 


Designation 

Date 

Origin 

Time 

N 

Lat. 

E 

Long. 

Depth 

mb 

Source 

NS  1 - WR  J3-23 

01-25-74 

12:27:47.  0 

58.  8 

27.  6 

- 

3.  0 

N 

NS1- WRUS-24 

01-27-74 

09:24:41.  0 

67.  7 

34.  3 

- 

3.  3 

N 

NS1  - WRUS-25 

01-27-74 

21:06:44.  0 

37.  5 

- 

3.  3 

N 

NS  1- FIND- 2 6 

01-29-74 

12:10:36.  0 

60.  5 

- 

3.  3 

N 

PS1  -GREC-27 

01-29-74 

15:12:45.  0 

38.  4 

31 

4.  3 

P 

j NS1-WRUS-29 

01-30-74 

14:43:20.  0 

- 

3.4 

N 

NS1- WRUS-30 

C 1-30-74 

14:50:34.  0 

59.4 

24.  8 

- 

3.5 

N 

PS1  -TURK-  31 

02-01-74 

00:10:62.  0 

38.  6 

27.  0 

29 

5.2 

P 

NS  1 - WR  US  - 3 2 

02-01-74 

14:06:17. 0 

59.  7 

27.4 

- 

3.  0 

N 

PS1  -TURK-33 

02-04-74 

23:38:04.  0 

27.  0 

57 

3.  6 

P 

NS1- WRUS-34 

02-06-74 

12:20:29.  0 

58.  6 

27.  9 

- 

3.2 

N 

PS1  -TURK-35 

02-07-74 

08:46:45.  0 

4 0.  1 

26.  9 

37 

4.  3 

P 

NS1  - WRUS-36 

02-08-74 

10:39:29.  0 

60.  3 

29.7 

- 

3.  3 

N 

NS1-WRUS-38 

02-13-74 

12:12:27.  0 

59.  9 

28.  7 

- 

3.  1 

N 

NS1  - VVRUS-4 0 

02-14-74 

13:41:18.  0 

60.  7 

29.  2 

- 

3.2 

N 

NS1  - WRUS-42 

02-15-74 

12:56:03.  0 

60.  0 

29.  2 

- 

3.  6 

N 

NS1-  WRUS-4  3 

02-15-74 

13:26:14.  0 

58.  2 

27.  5 

- 

3.2 

N 

NS1  - VVRUS-44 

02-15-74 

14:50:31. 0 

60.  8 

29.  3 

- 

3.  2 

N 

PS1  -GREC-46 

02-17-74 

05:03:07.  0 

38.  0 

21.7 

33 

3.  9 

P 

PS1  -IONI-4  7 

02-17-74 

11:33:15.  0 

37.  3 

20.  8 

8 

4.  0 

P 
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TABLE  II- 1 


EVENT  LIST  AND  PARAMETERS 
(PAGE  3 OF  3) 


Designation 

Date 

N 

Lat. 

E 

Long. 

Depth 

mb 

Source 

NS1-  WRUS-48 

02-18-74 

13:09:59.  0 

58.  7 

24.  0 

3.2 

N 

NS1-WR  US-50 

02-19-74 

12:38:53.  0 

58.  3 

27.  9 

- 

3.  5 

N 

NS1  -FIND-5  1 

02-21-74 

12:08:46.  0 

61.  3 

29.2 

- 

3.  3 

N 

NS1  - VVRUS-52 

02-22-74 

12:18:27.  0 

60.  8 

29.  3 

_ 

3.  5 

N 

PS1  -GREC-53 

02-23-74 

01:28:46.  0 

38.  1 

21.  8 

44 

4.4 

P 

PS1-GK  3C-54 

02-23-74 

20:55:19.  0 

38.  3 

20.  2 

33 

4.  4 

P 

PS1  -IONI-57 

02-26-74 

11:13:21.  0 

37.  2 

20.  8 

33 

4.4 

P 

PS1-NSVA-59 

02-28-74 

22:18:40.  0 

81.  6 

-2.  0 

33 

4.  5 

P 

PS1  -SVAL-60 

03-03-74 

16:22:48.  0 

76.  5 

13.  1 

33 

4.2 

P 

PS1-SVAL-61 

03-06-74 

00:10:00.  0 

77.  0 

12.  7 

33 

3.  9 

P 

PS1  -GREC-62 

03-10-74 

21:51:03.  0 

40.  8 

21.2 

17 

4.  1 

P 

PSI  - ALBN-63 

03-14-74 

20:57:34.  0 

41.  9 

19.4 

33 

4.  3 

P 

PS1-GREC-65 

03-22-74 

17:02:21,  0 

40.  7 

20.  5 

33 

4.  7 

P 

PSI  -S1CI-66 

03-28-74 

21:32:35.  0 

37.  1 

15.  0 

33 

4.6 

P 

PSI  - ALGE-67 

04-07-74 

00:56:56.  0 

36.  5 

04.4 

33 

4.  2 

P 

PSI  -ALBN-68 

04-08-74 

00:09:42.  0 

40.  6 

19.8 

19 

3.2 

P 

PSI  -GREC-69 

04-09-74 

00:31:09.  0 

39.4 

2 0.2 

44 

3.  1 

P 

PSI  -RUMA-7 1 

04-17-74 

01:31:34.  0 

46.  0 

21.  1 

33 

5.  6 

P 

PSI  -NNOR-73 

04-28-74 

12:52:50.  0 

68.  7 

16.2 

33 

4.  9 

P 

NS1 -URAL-74 

04-30-74 

14:17:15.  0 

57.7 

56.  9 

- 

3.6 

N 

P = NOAA-PDE  List  N = NORSAR  List  - = No  Information 
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Figure  II- 1 shows  location  of  NORSAR  and  the  event  epicenters,  indicated  by 
crosses.  The  circles  about  NORSAR  are  at  5 degree  increments.  Further, 
the  data  reported  by  Laun  and  Becker  (1974)  are  also  included  in  this  report. 

The  event  listing  for  their  data  is  given  in  Appendix  A. 

The  accuracy  of  the  event  parameter  information  (Table  II - 1 ) 
varies  considerably  between  source  bulletins.  The  NOAA-PDE  information 
is  generally  quite  reliable  since  many  stations  were  used  to  determine  epi- 
central  parameters  and  m^  On  the  other  hand,  for  the  NORSAR  bulletin,  the 
event  parameters  can  have  large  errors  due  to  beamforming  problems  associat- 
ed with  regional  events  (Ringdal  and  Whitelaw,  1973).  .Locations  errors  can 
be  as  much  as  160  kilometers.  Some  effects  of  these  problems  are  discussed 
in  the  following  sections. 
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+ Event  Locations 


FIGURE  II- 1 

STATION  AND  EVENT  LOCATION  MAP 


SECTION  III 


PROCESSING  AND  ANALYSIS 


A.  PROCESSING 

The  raw  NORSAR  short-period  data  were  edited  from  the 
library  tapes  from  Norway.  The  edit  proc  edure  reformats  the  single  sensor 
time  series,  performs  simple  quality  checks,  and  transfers  the  data  to  a 
working  tape.  This  edit  program  can  handle  192  seconds  of  data  at  a rate 
of  10  samples  per  second.  Two  edits  of  192  seconds  were  run  for  each  event. 
The  first  starts  about  20  seconds  before  the  P arrival  while  the  second  edit 
joins  the  first  to  form  384  seconc.s  of  contiguous  data.  Even  with  384  seconds 
of  data,  some  of  the  later  phases  for  the  more  distant  events  were  not  in- 
cluded in  the  edit  time.  The  more  distant  events  were  included  in  the  ensemble 
to  determine  the  distance  at  which  the  first  arrival  is  the  teleseismic  P and 

to  compare  regional  versus  teleseismic  m estimates. 

b 

Laun  and  Becker  (1974)  determined  the  sensor  which  had  the 
best  signal -to -noise  ratio  (SNR)  to  be  the  center  sensor  at  subarray  four. 

For  this  study  all  event  measurements  and  analysis  utilize  the  data  from  this 
seismometer.  From  inspection  of  the  noise  characteristics,  they  also  found 
that  the  raw  data  was  generally  dominated  by  two-second  period  noise  so  a high- 
pass  filter  (from  0.  9 to  5.  0 Hz)  was  applied  to  the  data.  A second  filter  was 
also  applied  to  correct  for  instrument  response.  Figure  III  - 1 is  a plot  of  the 
instrument  response  for  the  short-period  system  at  NORSAR. 

Phase  energy  ratio  determinations  were  made  using  the  Seismo- 
print  technique  (Cohen,  1969).  This  technique  displays  the  power  of  a seismic 
signal  as  a function  of  time  and  frequency  by  estimating  successive  power 
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spectra  from  a set  of  overlapping  time  gates  of  the  data.  The  spectral  power 
in  each  gate  is  normalized  to  the  largest  power  value  of  the  entire  set  and 
the  resulting  spectra  are  then  displayed  using  t wo -dimensional  printer  plots 
with  power  coded  by  letters.  All  spectral  estimates  used  a 32  point  time  gate 
(3.2  seconds),  contiguous  spectral  estimates  are  separated  by  8 points  (0.8 
seconds).  This  choice  was  considered  optimum  by  Cohen,  19f-9.  In  addition, 
the  program  computed  the  total  power  in  each  time  gate  for  use  in  the  phase 
energy  determinations. 

B.  VISUAL  ANALYSIS 

For  this  analysis,  a plot  of  all  subarray  sensors  was  examined 
to  check  for  continuity  of  the  chosen  P phase  across  the  subarray.  This  pre- 
vented picking  spurious  phases  at  the  reference  sensor.  In  most  cases  the 
P wave  was  clearly  visible  on  these  seismograms.  For  these  data  arrival 
times,  amplitudes,  and  periods  were  measured  for  estimation  of  location  ac- 
curacy and  determination  of  bodywave  magnitudes  m^. 
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SECTION  IV 
RESULTS 

A.  TRAVEL-TIMES  OF  P OR  P 

Ti  blc  IV-1  lists  the  observed  P travel-times  and  Figure  IV- 1 
shows  the  reduced  travel-times  relative  to  the  Herrin  (1968)  predicted  travel- 
times  for  zero  and  40  kilometers  depth. 

The  presumed  explosion  travel-times  scatter  about  the  predict- 
ed curve  wiih  deviation  of  about  2-3  seconds.  On  the  other  hard,  the  times 
for  earthquakes  located  by  the  PDE  are  greater  ‘han  those  predicted  for  nor- 
mal depth.  This  is  especially  evident  for  the  four  earthquakes  located  by 
PDE  at  distances  less  than  1800  kilometers  which  lie  near  the  0 km  depth 
predicted  times.  The  problem  presumably  lies  with  the  crust  and  upper  man- 
tle structure  (Archambeau,  et  al.  , 1969),  which  probably  does  not  conform  to 
the  average  world-wide  structure  inferred  from  the  observations  by  Herrin. 
The  earthquakes  located  by  NCRSAR  alone  scatter  significantly  about  the  Her- 
rin curve,  indicating  the  poor  location  capabilities  of  NORSAJt  (or  for  that 
matter,  any  single  array  station). 


B.  BODYWAVE  MAGNITUDE  (m  ) ESTIMATIONS  AT  REGIONAL 
’ DISTANCES 

J Bodywave  magnitude  is  an  important  parameter  in  many  tele- 

seismic  discriminants.  This  section  deals  with  the  relationship  between  m^'s 
| estimated  at  a single  NORSAR  instrument  at  regional  distances  and  those  esti- 

mated teleieismically  by  the  PDE.  The  fact  that  differences  exist  between 

I 

I 
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TABLE  IV  - 1 


OBSERVED  Pp  OR  P TR  A VFL  - TIMES  AND  AMPLITUDE) 

(PAGE  1 OF  3) 


Evcit 

■ 

1 

mSSm 

Phase 

T ravel  - 
Time 
Observed 
(sec) 

Period 
(sec ) 

All' 

(mil/  sec) 
Peak-to  - 
Peak 

mb 

NS1-WRUS-01 

r 

8.  88 

— i 

85.  9 

Pn 

125.  3 

0.  30 

12.  60 

4.  34 

PS1  - TURK-0 2 

23.  27 

147.  8 

P 

304.  4 

0.  95 

10  00 

3.  88 

NSl-U’RUS-0  3 

18.  22 

90.  7 

P 

248.  6 

0.  30 

10.  70 

3.  98 

PS1  -CITA-04 

17.  52 

176.  9 

P 

233.  4 

0.85 

3.  72 

3.  61 

PS1  -AEGE-0  5 

22.  33 

151.8 

P 

293.  2 

0.  85 

5.  75 

3.  t>4 

NSl-U'RUS-07 

8.  48 

90.  7 

P 

n 

120.0 

0.  60 

10.  30 

4.  13 

PS1  -AEGE-08 

22.  40 

152.  1 

P 

295.  3 

0.  90 

21.  90 

4.  23 

NS1  -SWRS-09 

20.  90 

1 15.  0 

P 

279.  5 

0.  90 

7.  57 

3.  68 

NS1-ADRI-10 

18.  26 

165.  9 

P 

250.  5 

0.  65 

7.  41 

3.  83 

PS1-SGRE-1I 

24.  33 

160.  4 

P 

315.  3 

0.  70 

5.  88 

3.  84 

NS1-WRUS-12 

8.  85 

86.  6 

Pn 

124.  8 

0. 90 

7.  42 

4.  01 

PS1 - TURK-  1 3 

23.  26 

143.  6 

P 

307.  4 

0.70 

8.  12 

3.  88 

NS1-WRIJS-14 

8.  39 

100.  3 

Pn 

119.  4 

0.  75 

9.  53 

4.  10 

NS1-WRUS-15 

8.  76 

92.  2 

Pn 

124.  7 

9.  85 

8.  91 

4.  08 

NS1-WRUS-16 

8.  89 

92.  6 

Pn 

119.  8 

0.  70 

10.  70 

4.  17 

NS1-WRUS-18 

8.  47 

92.  8 

Pn 

120.  4 

0.  oO 

7.  75 

4.  01 

PS1  -GREC-19 

23.  29 

161.9 

P 

306.  0 

0.  95 

91.  10 

4.  93 

NS1  - WRUS-20 

8.  97 

87.  1 

Pn 

127.  4 

0.  60 

12.  0( 

4.  23 

NS1  - WRUS-21 

8.  70 

95.  1 

Pn 

122.  2 

0.  40 

5.  47 

3.  86 

NS1- WRUS-22 

8.  87 

82.  0 

Pn 

125.  3 

0. 60 

2.  04 

3.  45 
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OBSERVED  Pn  OR  P TRA VFL-TIMFS  AND  AMPLITUDES 

(PAGE  2 OF  I) 


Distant  e 
Degrees) 


NSI-WRT"'-23  8.57 
NS1-WRUS-24  12.06 
NS1-WRUS-25  24.50 
NS1  - FIND -26  8.74 

PS1  -GREC-27  23.50 
NS1-WRUS-29  6.70 
NS1  - WRUS-30  7.00 
PS1-TUKY-31  24.40 
NS1-WRUS-32  8.20 


PS1  -TURK-33  22  90  147.  3 


NS1-VRUS-34  8.80 


NS1  - WRUS-36  9.20 

NS1  - WRUS-38  8.80 


NS1-WRUS- 


8.  90 


NS1  - WRUS-42  9.00 


NS1  - WRUS-44  8.90 


97. 7 


PS1  -TURK-35  23.00  147.6 


85.  5 
88.  8 
83.  3 
87.  8 


NS1-WRUS-43  8.80  100.7 


82.  6 


PS1  -GREC-46  23.90  158.8 

PS1  -IONI-47  24.40  160.9 


0.  70 


318.  8 0.  85 


A / T 

( m ji ' set  ) 
Peak-to- 
Peak 


3.  89 

4.  28 

3.  95 

4.  02 
4.  23 


3.  92 

3.  57 

4.  17 


4.  22 
4.  17 
4.  06 
4.  27 


7.23  3.99 


14.69  4.20 


r 


TABLE  IV-1 

OBSFRVED  Pn  OR  P TRAVEL-TIMES  AND  AMPLITUDES 
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Event 

Di  stance 
(Degrees) 

A zi  muth 
STA- 
EPI 

Phase 

T ravel  - 
Ti  me 
Observed 
(set ) 

Period 

(sec) 

A/T 

(m/i/set  ) 
Peak  -to  - 
Peak 

m 

b 

NS1  - WRUS-48 

6.  90 

102.  9 

Pn 

- 

0.  30 

3.  74 

3.  57 

NS1-WRUS-50 

8.  90 

99.  5 

Pn 

129.7 

0.  40 

12.  90 

4.  25 

NS  1- FIND -51 

8.  80 

79.  4 

Pn 

127.  9 

0.  50 

4.  67 

3.  81 

NS1-WRUS-52 

8.  90 

82.  6 

Pn 

129.  6 

0.  50 

12.  31 

4.  23 

PS1  -GREC-53 

23.  80 

158.  6 

P 

309.  3 

0.  95 

18.  60 

4.  29 

PS1-GREC-54 

23.  30 

161.  5 

P 

306.  1 

0.  85 

10.  70 

4.  00 

PS1-IONI-57 

24.  50 

160. 9 

P 

317.  3 

0.  85 

7.  08 

3.  94 

PS1-NSVA-59 

21.  00 

354.  7 

P 

282.  7 

0.  85 

38.  00 

4.  38 

PS1-SVAL-60 

15.  60 

1.  8 

P 

218.  4 

0.  40 

13.  49 

4.  36 

PS1-SVAL-61 

19.  70 

1. 4 

P 

226.  2 

0.  90 

6.  76 

3.  61 

PS1-GREC-62 

21.  10 

158.  1 

P 

283.  7 

0.  85 

9.  76 

3.  79 

PS1  -ALBA-63 

19.  70 

161. 2 

P 

- 

0.  85 

3.  89 

3.  86 

PS1  -GREC-65 

21. 00 

159.  6 

P 

283.  5 

0.  80 

5.  62 

3.  55 

PS1  -SICI-66 

23.  90 

172.  1 

P 

313.  2 

0.  70 

3.  99 

3.  62 

PSl  -ALGE-67 

24.  80 

192.  7 

P 

320.  5 

0.  90 

7.  41 

3 99 

PS1-AI  BN -68 

21.  00 

161.  1 

P 

284.  0 

0.  50 

10.  50 

3.  82 

PSl  - GREG -69 

22.  20 

161. 0 

P 

294.  5 

0.  Q0 

5.  24 

3.  59 

PC1  - RUM A - 71 

16.  00 

153.  8 

P 

225.  8 

0.  80 

3.  46 

3.  74 

PS1-NNOR-73 

8.  10 

13.  5 

Pn 

118.  4 

- 

- 

- 

NS1  - URAL  -74 

23.  20 

77.  4 

P 

310.  5 

0.  60 

5.  37 

3.  78 

REDUCED  P TRAVEL-TIMES 
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regional  estimates  and  tho  je  for  teleseismic  distances  is  well  known 

(Bungum  and  Husebye,  1974).  For  this  study  NOAA-PDE  values  computed 

utilizing  stations  at  distances  less  than  20  degrees  were  routinely  recomputed 

omitting  the  close-in  stations  in  order  to  ensure  that  they  were  the  'true'  tele- 

seismic  m,  's. 

b 

Figure  IV-2  shows  the  NORSAR  single  instrument  m estimates 

b 

minus  the  reported  NOAA-PDE  or  NORSAR  beamed  m 's.  We  show  these  data 

L' 

as  a function  of  distance  to  ascertain  whether  the  standard  distance  factors  (B°) 

used  were  suitable  for  the  calculation  of  m,  in  Eurasia.  However,  the  scatter 

b 

cf  the  data  is  too  large  to  make  any  conclusion  concerning  this  problem. 

It  is  apparent  that  there  are  two  populations  formed  by  this  sin- 
gle subtraction.  The  NORSAR  single  instrument  m^'s  are  larger  than  those 
determined  from  the  NORSAR  beams  by  about  0.8  magnitude  units,  while  those 
determined  from  NOAA-PT  £ are  larger  than  the  single  instrument  NORSAR 
m^'s  by  about  0.  3 to  0.  4 magnitude  units. 

Bungum  and  Husebye  (1974)  indicate  a bias  of  about  0.2  magni- 
tude units  between  the  teleseismic  NORSAR  beam  and  NOAA  PDE  m,  's.  This 

b 

is  compatible  with  the  bias  found  here.  Further,  they  found  that  the  bias  for 
the  NORSAR  beamed  m^'s  for  various  distances  and  regions  are  dependent  on 
event  magnitude.  These  results  are  mostly  due  to  the  threshold  bias  problem 
which  arises  when  PDE  magnitudes  are  calculated  by  averaging  over  detect- 
ing stations  only  (Ringdal,  1975  ).  A plot  of  NOAA-PDE  magnitudes  versus 
NORSAR  single  site  magnitudes  confirms  this  threshold  bias  while  suggesting 
that  some  bias  may  be  due  to  propagation  effects.  The  data  sample  is  too 
small  to  verify  the  latter  conjecture. 

The  NORSAR  beamed  estimates  are  comprised  of  many  sm,i!l 
magnitude  events  (m^S  3.6)  and  the  large  NORSAR  beam-single  sensor  dif- 
ference (0.8  magnitude  units)  is  due  primarily  to  beamforming  losses  inten- 
sified by  the  presence  oi  relatively  high  frequency  energy.  This  high  frequency 
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AS  A FUNCTION  OF  DISTANCE 


energy  increases  beamforming  losses  because  of  greater  signal  dissimilarity 
at  high  frequencies  across  the  array. 

In  Figure  IV-3  the  high  frequency  content  is  clearly  visible  in 
trace  1,  a single  sensor  high-pass  filtered  NORSAR  trace,  just  after  the  sig- 
nal arrival  time.  In  fact  this  signal  would  not  have  been  detected,  were  it  not 
for  the  presence  of  this  energy.  Trace  2 shows  the  same  time  series  correct- 
ed for  the  NORSAR  instrument  response  and  no  signal  is  present.  Hence  an 
instrument  peaking  around  1.  0 Hz  would  probably  no*  have  detected  this  event. 
The  NORSAR  instrumen’  pu  iks  ai  3.  0 Hz,  and  is  ideal  for  the  detection  of 
small  regional  events. 

Summarizing,  we  find  a negative  bias  of  0.2  to  0.4  magnitude 
units  for  and  slow  tr  avel-times  for  regional  events  observed  at  NORSAR. 
This  suggests  low  velocity  and  possibly  high  attenuation  in  the  crust  and  up- 
per mantle  along  the  propagation  paths  to  NORSAR,  although  this  result  is 
only  suggested  by  the  data. 

The  scatter  of  single -instrument  NORSAR  estimates  pre- 
vents the  determination  of  a conventional  constant  station  correction.  Sec- 
ondly, the  relationship  between  the  NORSAR  beamed  m^  estimates  for  the 
small  events  (m^£3.6)  and  the  teleseismic  magnitudes  is  unknown.  This  re- 
lationship can  not  be  extrapolated  from  these  data  and  can  not  be  determined 
experimentally,  due  to  the  lack  of  corresponding  PDE  magnitudes. 

C.  P/S  ENERGY  RATIOS 

As  discussed  earlier,  many  of  the  locations  were  not  well  de- 
fined and  hence  individual  phases  were  often  not  identified  due  to  timing  errors. 
That  being  the  case  comparison  of  amplitudes  and  times  between  phases  might 
be  misleading  and  we  concentrate  on  gross  differences  between  the  two  sources 
as  did  Booker  and  Mitrinovas  (1964). 
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The  short-period  seismogram  was  divided  into  three  broad 
time  frames  to  take  advantage  of  the  differences  between  earthquakes  and  ex- 
plosions. These  three  time  frames  contain  predominantly  compressional  wave 
energy,  shear  wave  energy,  and  surface  wave  energy,  in  that  order,  and  are 
based  on  approximate  measured  travel-time  using  the  following  velocity  win- 
dows : 

• P 8,  6 km/ sec  to  8.  0 km/  sec 

• S 5.  0 km/ sec  to  4.  Okm/ sec 

• Rg  4.0  km/sec  to  3.  Oku/sec. 

The  total  power  in  each  spectrum  computed  by  the  seismoprint  program  was 
summed  over  the  times  corresponding  to  these  velocity  ranges,  and  the  ratios 
between  them  computed.  Because  the  short-period  record  length  was  limited 
to  384  seconds,  the  complete  Rg  phase  for  A > 15°  was  not  available  and  the 
complete  S phase  for  A > 20°  was  also  not  available. 

Figure  IV-4  shows  the  P/S  energy  ratios  as  a function  of  dis- 
tance. Several  general  observations  can  be  made: 

• There  is  separation  using  this  criterion  between  presumed  ex- 
plosions and  earthquakes  but  there  are  too  few  presumed  explo- 
sions to  state  any  conclusions  with  confidence.  Laun  and  Becker 
(1974)  showed  the  same  presumed  explosion  data  and  it  is  unfor- 
tunate that  no  additional  presumed  explosions  were  available  for 
analysis  during  the  timr‘  span  of  our  event  ensemble. 

• P/S  ratios  for  earthquakes  having  epicenters  determined  by  PDE 
are  more  consistent  (i.  e.  , have  less  scatter)  than  for  those  lo- 
cated by  NORSAR.  However,  there  are  only  5 such  values  which 
are  certainly  not  enough  for  a good  statistical  sample. 
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• Laun  and  Becker  discussed  a presumed  explosion  and  earthquake 
which  occurred  within  one  degree  of  each  other  in  western  Rus- 
sia and  about  1700-1800  kilometers  distance  from  NORSAR.  Our 
ensemble  of  data  added  one  earthquake  (PS1  -RUMA-7  1 ) located 
in  Rumania,  whose  P/S  ratio  was  similar  to  that  of  the  western 
Russia  earthquake. 

• The  presumed  explosions  occurring  at  distances  of  about  2180 
kilometers  arc  from  Novaya  Zemlya.  All  of  them  have  about 
the  same  P/S  ratios.  However,  there  are  insufficient  earth- 
quake data  to  make  any  statements  concerning  the  degree  of  sep- 
aration between  source  types  at  this  distance. 

It  is  clear  that  we  need  more  presumed  explosion  data  (P/S 
ratios)  to  determine  with  conficence  what  the  degree  of  separation,  if  anv,  is 
possible  with  this  criterion. 
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SECTION  V 


SUMMARY 


In  the  basic  analysis  of  small  regional  events  observed  at 
NORSAR,  we  have  looked  at  P travel-times,  bodywave  magnitudes  and  one 
promising  discriminant,  short-period  vertical  P/S  energy  ratios. 

Errors  in  epicenter  locations  resulted  in  erroneous  station  to 
epicenter  distances  and  consequently  analysis  of  phases  other  than  P were 
questionable.  mfa  estimates  are  biased  by  the  presence  of  relatively  high  fre- 
quency energy,  erroneous  distances,  and  upper  mantle  structural  effects. 

P/S  energy  ratios  appear  promising  as  a discriminant  even 
with  this  poor  experimental  setup.  Ideally  one  should  have  a network  of  sta- 
tions observing  seismic  events  at  regional  to  small  teleseismic  distances. 

With  such  an  experimental  plan,  event  locations  (hypocenters)  could  be  made 
with  sufficient  accuracy  to  conduct  a complete  first-zone  study.  Unfortunately, 
such  data  for  Eurasia  was  not  currently  available. 
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APPENDIX  A 


EVENT  LIST  AND  EVENT  MEASUREMENTS 
FROM  LAUN  AND  DECKER  (1974) 
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EVENT  LIST  AND  PARAMETERS 


Abbreviations: 


MAGNITUDE  MEASUREMENTS  FROM  THREE  SOURCES 


LAUN  AND  BECKER  (1974) 


Event 

Number 

NOAA-PDE 

m, 

b 

NORSAR 

Beam  in, 
b 

NORSAR 

Single-Site  m, 
b 

1 

4.  5 

ND 

4.  31 

2 

4.  6 

ND 

4.  43 

3 

- 

3.  1 

4.  05 

4 

- 

3.  2 

4.  21 

5 

- 

3.  5 

3.  88 

6 

- 

3.  4 

3.  90 

7 

- 

3.  3 

3.  93 

8 

- 

3.  5 

3.  75 

9 * 

6.  8 

ND 

6.  50 

10  * 

6.  0 

ND 

5.  59 

1 1 

- 

3.  5 

3.  90 

12  * 

5.  2 

ND 

5.  19 

1 3 

- 

3.  3 

ND 

14 

- 

3.  3 

3.  78 

15 

- 

3.  2 

4.  11 

16  * 

4.  8 

4.  9 

5.  29 

17  * 

6.  9 

ND 

5.  45 

18  * 

4.  2 

3.4 

3.  46 

19  * 

4.  4 

3.  1 

3.  40 

20  * 

4.  8 

3.  8 

3.  82 

- - . . „ 

■t*  Recomputed  using  teleseismic  m values  only 

b 

* Presumed  explosion 
ND  Not  Detected 
No  value  given 
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